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Abstract 



An interesting feature of spin- 1 chains with a gap is that they undergo a 
commensurate-incommensurate transition in the presence of an external mag- 
netic field H. The system is in a gapless incommensurate regime for all values 
of the magnetic field between the lower critical field H c \ and an upper critical 
field H C 2, where it is gapless and has power law correlations. We calculate the 
critical exponents for such a generic gapped system in the incommensurate 
regime at the critical field H c \ and in its vicinity. Our analysis also applies 
to the spin-^ ladder. We compute the full dynamical susceptibilities at finite 
temperature. We use the same to discuss the thermal broadening of various 
modes and obtain the low temperature behaviour of the nuclear spin relax- 
ation rate T£~ . We discuss the results obtained here for the special cases of 
the dimerised chain, frustrated chain and the spin-^ ladder. 



I. INTRODUCTION 

Quantum spin chains have been a field of intense theoretical activity for the past decadeS. 
This seems pertinent in the face of the tremendous advance made in the fabrication of 
quasi-one dimensional systems. There are a lot of compounds like (TMTTF) 2 X which are 
essentially spin-| chains, NENP - a spin-1 chain! and a host of other organic compounds 
whose behaviours can be adequately described within the framework of interacting spin 
systems!. The Heisenberg model with purely nearest neighbour interactions has been in- 
strumental in understanding various properties of these spin systems@. In addition to these, 
there are compounds with alternation in the nearest neighbour bonds (arises as a result of 
interactions with phonons) eg TMTTF, strong single ion anisotropies etc. The effects of 
spin-Peierls dimerization, interchain interactions and competing nearest neighbour interac- 
tions have been reasonably well understood in the spin-| caseO. These are known to drive 
phase transitions in the pure spin-| model. The basic issue addressed is whether there exists 
a gap above the singlet ground state to the first excited state which is a triplet. Apart from 
these, the ground state could also exhibit crossovers from short range Neel order to spiral 
order or spontaneously dimeriseQ but these are not signals of any kind of quantum phase 
transitions. 

More interesting phases can be obtained by placing these gapped systems in a magnetic 
field. For some critical value of the field H&, the system undergoes a continuous phase 
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transition from a commensurate Neel (C) zero uniform magnetisation phase to an incom- 
mensurate phase (IC) with non-zero magnetisation. As the magnetic field is increased, the 
magnetisation increases and saturates at some critical H c2 where the ground state is fully 
ferromagnetic but differs from the XXX ferromagnet in that it has a finite gap to the first 
excited state. The intermediate region (between H c i and H c2 ) is completely gapless and 
the pitch vector Q (the value of the momenta at which the static structure factor shows 
a peak ) decreases continuously from Q = n at H c \ to Q = at H c2 . Evidence for such 
transitions were seen in the spin-| ladder compound Cu 2 (C 5 H 12 N 2 ) 2 Cl$ and the quasi one- 
dimensional spin-i systems TTF — CuBDTu. The C-IC transition seen in 
can be understood in a similar manner if it is treated as a spin chain with spin-Peierls inter- 
actions alone. However, because of the presence of strong phonon interactions in CuGeO^, 
various complications can arise, in particular, the transition becomes first-order. 

In this paper, we are concerned with the properties of a generic gapped spin-| system 
in the presence of an external magnetic field. Though the naive expectation is that all 
gapped systems might exhibit similar behaviours in the presence of the field, we find that 
they have very different properties depending on the nature of the interaction which creates 
the gap. In particular we focus on experimentally measurable quantities like the total 
magnetisation, NMR relaxation rates and neutron scattering intensities. To do this we 
calculate various spin-spin dynamic correlation functions in the gapless region between H c \ 
and H c2 . We calculate various exponents in the IC regime close to H c \ and use these to 
calculate the temperature dependence of the NMR rates (Ti). We also discuss in detail 
the consequences of our results for the dimerised, frustrated and ladder systems. We find 
that though magnetisation measurements will not help differentiate between these systems, 
neutron scattering and NMR can. Since the models studied in this paper directly describe 
the compounds mentioned above, the results obtained here are of immediate relevance to 
experiments. 

The paper is organised as follows: Sec. II contains a very brief review of spin-^ chains and 
the technique of bosonization used in this paper. In Sec. Ill , we introduce our hamiltonian 
for a generic gapped spin-| chain and discuss the effects of the applied magnetic field and the 
consequences of incommensurability induced by the magnetic field. We present in Sec. IV 
our results for the dynamic correlations and the values of the exponents near the transition 
at H cl for a generic gapped spin-| chain. We also discuss the generic phase diagram as a 
function of magnetic field and temperature. The spin ladder is treated separately in Sec. 
V. In Sec. VI, we calculate the temperature dependence of the NMR rate T{~ and discuss 
the behaviour of the same for the cases of the dimerised, frustrated and ladder models. Sec. 
VII contains the concluding remarks. 



II. INTRODUCTION TO SPIN-± SYSTEMS 

Here we present a brief introduction to the physics of spin-| chains and summarise 
certain results relevant to the work presented here. The main hamiltonian of interest is the 
Heisenberg model: 

H = JY,Si-S i+1 (1) 
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where the Si is a localised spin | operator. We set J = 1 in the rest of the calculation. This 
model is Bethe ansatz soluble and is known to be gapless. A lot is known about the static 
and dynamic properties of this system. Eventhough there are various methods used to study 
spin-| chains, in this paper we use the machinery of bosonization since this enables us to 
deduce various properties in a relatively easy manner. Since this is an oft used method we 
sketch the details briefly. We first use the Jordan- Wigner transformation&Eil which essentially 
maps the spin problem onto a problem of interacting fermions on a lattice. For the spin-| 
system considered here, the corresponding fermionic problem has fermi momentum &f = |. 
We then perform a linearisation around the free fermi points given by ±kp , to obtain 
an effective low energy continuum fermionic theory and then bosonize using the standard 
dictionary of abelian bosonizationEH. We just present the final expressions obtained for the 
spin | operators in terms of the bose field <fi and its dual 9 

^(x) = (-ircos(20) + -U x 

S~(x) = exp{i6)[{-l) x + cos(20)] (2) 

Finally, Eqs.(2) can be used to obtain the bosonized version of the hamiltonian given in 
Eq.(l) i.e., 

H = h I dx ^ uK ^ 2 + (^)(^) 2 ] ( 3 ) 

where u is the spin wave velocity, IT is the momentum conjugate to the field and K = | 
for the isotropic hamiltonian of Eq.(l). 7i is just the hamiltonian for free bosons. Note that 
K = 1 for the case of the XX- antiferromagnet which in turn is equivalent to a theory of free 
fermions via the Jordan- Wigner transformation. Other values of K correspond to having a 
J z coupling and hence to interacting fermions. Since there is no mass term for the field, it 
is clear that there is no gap to the first excited state from the singlet ground state. Since a 
free boson theory given by Eq. (3) is trivially solvable, it is fairly straightforward to calculate 
the dynamic correlation functions using Eqs.(2). 

III. SPIN-i SYSTEMS WITH A GAP IN A MAGNETIC FIELD 

A. Zero Magnetic Field Case 

In many compounds which mimic the behaviour of spin-| chains, it has often been found 
that there are many other interactions between spins, apart from the the isotropic inter- 
action summarised in the hamiltonian of Eq.(l). These could be anisotropies, frustration, 
impurities and interchain couplings. More often than not, it is found that these interactions 
are small and can be treated as perturbations of the above gapless system. Within the 
framework presented here, the effects of these perturbations can be gauged from whether 
they are relevant, marginal or irrelevant in the renormalisation group sense. Depending on 
whether the perturbations are relevant or irrelevant, a gap opens in the spectrum. In this 
paper we are interested in a weak spin-Peierls dimerisation and frustration arising from next 
nearest neighbour interactions. 



3 



We first consider a small dimerization S, which arises from spin-lattice interactions. The 
corresponding perturbation to the hamiltonian given in Eq.(l) is 



^ = (-1)^.^+1 (4) 
Using Eqs.(2) the bosonized continuum version of (4) is§ 

H D = 5 J dxcos(20) (5) 

It is equally interesting to study the effect of a competing next nearest neighbour interaction 
J 2 . Such an interaction tends to frustrate the system. The addition to Eq.(l) is 

Hf = J^Si ■ S i+2 (6) 

For classical spins, the ground state retains its Neel like order for all J 2 < 0.25 and exhibits 
spiral order for J 2 > 0.25. This is just a crossover as it is not possible for J 2 to drive a 
zero temperature phase transition in the classical system. What is the behaviour of the 
quantum spin-| system ? To study this we again take recourse to Eqs. (2) and write down 
the bosonized version of Eq. (6)§ 

U F = (J 2 - J 2c ) J dx cos(40) (7) 

Apart from the cosine terms mentioned above, these spin interactions also renormalise 
the velocity u in (3). We now analyse the effects of the interactions given by Eqs. (6,7). 
In the free boson theory of Eq.(3), the anomalous dimensions of the operators cos(n</>) are 
given by IL ^. Using this at the isotropic point K — |, we see that Hd given by Eq.(5) 
has dimension | and hence is relevant for all values of S. Similarly Tip described in Eq.(7) 
is found to have dimension 2 and is marginal. Numerically it has been shown that 7ip is 
marginally relevant for J 2 > J 2c = 0.2411 and irrelevant for J 2 < 0.24110. Therefore, we 
can see that both S and J 2 drive a quantum phase transition from a gapless phase to a phase 
where a gap opens in the dispersion for all values of 5 and J 2 > J 2c . 

Instead of restricting ourselves to the dimerised and frustrated models in this paper, we 
study the following hamiltonian generic to gapped spin-| chains in a magnetic field: 

If u 
H g ap = — J dx[uK(TiYi) 2 + (-)(d x <pf + vcos(n<P)} (8) 

The 2-parameter space spanned by K and n correspond to various spin-| systems with 
gaps in their spectra. For example, K — ^,n — 2 corresponds to the dimerised model and 
K = |,n = 4 is the frustrated antiferromagnet. 

We see that though all these models have a gap, we can differentiate between them 
by studying their excitation spectra the nature of which depends on the values of K and 
n. For instance, it is known from the study of excitations in the sine-Gordon theories^!, 
that apart from the magnons and solitons there exist certain excitations called breathers. 
These are non-linear like the solitons and manifest themselves as discrete levels in excitation 
spectrum. These breathers have energies higher than the singlet-triplet gap but lie below the 



4 



continuum part of the spectrum generated by the solitons. The number N of such breather 
modes is given by iV < (8 — n 2 K)/n 2 K. For example, we see that the frustrated model has 
no breather modes in its excitation spectrum. On the other hand the spin-Peierls system 
has two breather modes in its spectrum: the N = 1 breather is just a renormalized magnon 
and the N = 2 breather is a bound state of two magnons. This extra bound state is found 
to affect the spin correlation function at finite temperatures^. Optical experiments on the 
lines of Ref. [T5] where breathers in quasi-lD ferromagnets were detected, can be used to 
directly detect these breathers. 



B. Effects of a Magnetic Field 

We now study the effect of a magnetic field on these gapped systems. To do so, we turn 
on a magnetic field H in the z-direction. Note that the magnetic field breaks the SU(2) 
symmetry. The interaction with the spins on the lattice is 

7i m = , £gn B HSt (9) 

i 

A simple understanding of the effect of the field can be obtained in the fermionic picture. 
First of all, the application of the magnetic field is equivalent to a chemical potential for 
the fermions. When H = we have zero uniform magnetisation i.e., J2iSf — 0. This 
corresponds to a completely filled lower band for fermions with a gap A separating the 
lower and the upper bands and the fermi energy lying in the middle. As we increase H, we 
are in effect shifting the fermi energy. A point is reached where the fermi energy crosses 
the gap and lies at the bottom of the upper band. This value of H corresponds to H c i i.e., 
A = gfisHd. When H is increased beyond H c i, the upper band is partially filled resulting 
in a non zero magnetisation m. The zero magnetisation region with a gap below H c i is 
the commensurate (C) regime in that, short range antiferromagnetic order still persists. 
Above H c i, the ground state is magnetised, has no gap and canted. This is what we call 
the incommensurate (IC) phase. This is a quantum phase transition and S z is the relevant 
order parameter. This is what makes it different from the incommensurate (spiral) phases 
seen in frustrated systems (Sec. IIIA)llHHl. 

Now that we have a heuristic understanding of the effect of the magnetic field we now 
proceed to calculate various quantities of interest like the magnetisation and exponents in 
the vicinity of H c i, within the continuum approximation. Using Eqs.(2) we can see that 
Eq.(9) corresponds to adding a gradient term gfiBHd x <f) to the hamiltonian of Eq.(8). 



Titot = 7-Lga P + 7T dxg/j, B Hd x (j) (10) 

This hamiltonian has the same form as the one used in the context of the C-IC transition in 
two-dimensional systems where it is a transition from a phase with a a discrete symmetry (C) 
to a phase with a continuous symmetry (IC)El'lil. Note that the gradient term in Eq.(10) can 
be eliminated by a simple shift of the <fr field i.e., — > <p + ^Kg^sHx. Since the only effect 
of H below if ei is to renormalise the gap, we can replace this shift by <fi — > <fi + rcmx where 
m is the magnetisation. The cosine term, however, is not invariant under this shift. From 
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the analogy with the chemical potential in fermionic systems, we infer that the magnetic 
field changes the fermi momentum. We can redefine new fermi points ±k' F , linearise around 
these points and obtain a new effective massless free boson theory, albeit with a different 
value of K. It has gapless modes at q = and q = 2k' F — 7r(l — 2m). A schematic plot of 
the field dependent dispersion for the spins in the IC phase is given in Fig.l. The behaviour 
in the gapless IC region is governed by the hamiltonian 

H = ±- I dx[vK^Yif + [t){d x <pf] (11) 
In J K 

The quantities v and K are dependent on H. Although it is difficult to obtain the depen- 
dences for the entire range of H between H c \ and one can compute the exact values for 
the exponents (equivalently the K) at the critical point H c \ and see how they change as we 
move away from these points. 

Before we proceed with the calculation of K, we first study how the magnetisation rises 
above H c \. Using the results on the C-IC transition!!! 



m = r) J (H 2 — H ( 



n 



4 2 i) 

= V^\j{H - H cl ){H + H cl ) ~ V^\j{H - H cl )2A (12) 

where n is the coefficient of (f> in the argument of the cosine interaction and rj is a constant of 
proportionality which depends on the renormalised velocities. The magnetisation increases 
from its zero value as a square root near H c \. A similar square root behaviour is seen 
near H c2 if one approaches this critical point from the ferromagnetic side. Here it is the 
decrease in magnetisation from the full ferromagnetic value that shows the square root 
behaviour. In Fig.2, we have shown the behaviour of the magnetisation with temperature 
and the field in units of g\iB- This was done using the analogy with the case of fermions 
with a gap and a chemical potential. We assumed a dispersion for the fermions of the form 



Uk = ±y(Jk) 2 + (A) 2 . Here the two square root regimes near H cl and H c2 (0.4 and 1.1 
in the figure) are joined by a region in which the magnetisation increases in a nearly linear 
fashion with H. The interval of H in which the square root behaviour of the magnetisation 
is seen and also the range where the nearly linear regime is seen, depend on the parameters 
like the couplings and the gap. For certain ranges of these parameters where the gap A is of 
the order of the exchange couplings one obtains a scenario where the width of the square root 
regime is quite small so that even experiments done at reasonably low temperatures could 
entirely miss the detection of this region. This might be the reason as to why the square 
root behaviour is not observed in Ref. |f[ Since the behaviour of the magnetisation given by 
Eq.(12) is generic to all the gapped systems studied here, it follows that a measurement of 
the magnetisation will not be able to qualitatively differentiate between the various models. 

The magnetisation derived in Eq. (12) applies to the ladder also. This is because, as 
we will later show in Sec. V, the effective hamiltonian for the ladder has the same form 
as given in Eq.(10). Another system which exhibits a similar transition in the magnetic 
field is the spin-1 chain. Since the spin-| ladder and the spin-1 chain belong to the same 
universality class^, Eq.(12) describes the magnetisation in the spin-1 chain also. This agrees 
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well with the result for m in spin-1 chains obtained in Ref. |20], [21} In Ref. |21] a mapping 
of a phenomenological hamiltonian of interacting magnons onto a system of bosons with 
repulsive 5-fn interactions was used to obtain the magnetisation and the correlations in the 
IC regime. 

Naively one might expect that these gapped systems in a magnetic field have the same 
qualitative behaviour for correlation functions and that quantitative features like H Cl and 
thermodynamic quantities like m are model dependent. We shall see below that this is not 
true and that these models have dynamically different physical behaviours in the IC regime 
depending on the value of n. For a generic interaction of the form cos(n0) (where n 2 K < 8 
) and H close to H cl on the IC side,E3 

R = 4 nm 7 sinh(20) 
n 2 A 



Here A is the gap and 7 is a positive constant which depends on the parameters of the 
theory but is independent of H. 9 is defined by exp(— 29) = At the transition H c i, 

K goes to a universal value \ independent of the value of K. As H increases, there is a 
non-zero magnetisation m and the change in K is completely governed by n and K in that K 
could increase or decrease depending on the value of K. For example, for the frustrated ( J 2 ) 
model K — \ at H — H c i and increases as H increases and for the case of the dimerised (S) 
model, K = 1 at H = H c \ and decreases with increasing H. Here we have seen that though 
the magnetisation has the same qualitative behaviour for all the models, the value of K is 
model dependent because it is determined by the anomalous dimensions of the perturbing 
operators. As a consequence the exponents in these models will be radically different as will 
be shown in the following section. 



IV. CORRELATION FUNCTIONS IN THE CRITICAL REGION 

In this section, we first study how the structures of the correlation functions are altered 
by the magnetisation in the IC phase. Later we will calculate the various commensurate and 
incommensurate contributions to the finite temperature dynamic susceptibility. From the 
magnetic interaction given by Eq.(9) and the field shift <fi — > <p + ttitix, it is clear that since 
for H < H c i there is no net magnetisation, the form of the correlators is unaffected and only 
the gap is renormalised. However, in the IC phase, the presence of a non-zero magnetisation 
m results in <ft —>■ <fi + Tcmx. The dual field 9 is insensitive to this shift. Incorporating this 
field shift in the expressions for the spin operators given in Eq.(2), we find that the generic 
form of the spin-spin correlation function in the IC region is 

(S z (x,t)S z (0,0)) = m 2 + /1OM) + cosvr(l - 2m)f 2 (x,t) 
(S + (x, t)S~(0, 0)) = cos(27imx)gi(x,t) + cos(7rx)g2(x,t) (14) 

where fi, / 2 , gi and g 2 are monotonic, decreasing, power- law functions of space and time 
and they go to zero asymptotically. From the above expressions we deduce the following: 
the correlation function parallel to the field i.e., (S Z S Z ), has a uniform magnetisation or 
equivalently a Q = mode while the staggered part is shifted from Q = 7rtoQ = 7i — 2^m. 
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This is in contrast to the correlation in the plane perpendicular to the field, (S + S~), where 
the staggered mode is unshifted and remains at Q = 7r and the uniform magnetisation (Q = 
0) mode is shifted to Q = 2nm. We call Q = 2nm and Q = tt(1 — 2m) the incommensurate 
modes. 

With the results obtained above , we now proceed with our calculation of experimentally 
relevant quantities like the dynamic susceptibilities at finite temperatures. As a first step we 
compute the associated unequal time correlation functions using Eqs.(2),(14) and the result 
of Eq.(13). The correlation functions are found to be as follows: 

(S z {x,t)S z {0)) = ({S z {x,t) - m){S z {0) - m)) 

K 1 1 

= cos(vrx(l - 2m))(x 2 - t 2 )~ K + cste. — (- — + 



'4vr 2V (x-t) 2 (x + tf 
(S + (x,t)S-(0)) = (-l) x (x 2 -t 2 )-^ + cste. cos(2vrmx) {x 2 - t 2 )-^ +k ~ l) ■ 

[exp(2i7rmx)- — + exp(-2inmx)- --} (15) 

[x — ty (x + ty 

with K being specified by Eq.(13). It is not surprising that the exponents are different in 
the directions parallel and perpendicular to the external field. This is because the magnetic 
field breaks the SU(2) spin rotational invariance. For the case of the dimerised model at the 
critical magnetic field H = H&, we found in Sec. Ill that K = 1. An interesting coincidence 
is that the exponents and hence the correlation functions calculated here for the dimerised 
model at H = if cl are the same as that for the XX- antiferromagnet in a zero magnetic 
field. The correlation functions of the XX-model (which is in turn equivalent to a theory of 
free fermions) can be obtained by substituting m = and K — 1 in Eqs.(15). Note that 
this correspondence holds only at the critical point. 

Using the above results we can also evaluate the (g, to) dependent susceptibility at finite 
temperatures, both in the direction of the applied magnetic field and in the direction per- 
pendicular to it. A knowledge of this quantity helps us extract various measurable quantities 
like neutron scattering intensities, absorption and nuclear magnetic resonance (NMR) rates. 
The susceptibilities are given by the following expression 



Xij{q,uj,T) = -i J dt dx expi(iot - qx)6(t)([S l (x,t) } S J (0,0)}) T (16) 

Here i, j refer to the components of the spin and the subscript T implies that the correlator is 
evaluated at finite temperature. Since rotations in the x — y plane still leave the hamiltonian 
invariant, there are no cross correlations i.e., (S 1 *^- 7 ) = for i ^ j. For the susceptibility 
Xzz = X|| in the i-direction parallel to the direction of the applied 

n 2 

Q=Q, rp\ H 

For the perpendicular susceptibility x± = X+~ i-e., in the x — y plane 
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xf— < s ,«,r) = -n>tV[b(^ - i ( " + ^ g)) , (-1 - «)*(§ - i ( ""^- q)) , (1 

where 2/3 = 2K + ± - 2 

xfto,«.r) = M*-B(f - (i - „) )B (| - ^-^-o)) , (1 _ o)) 

(19) 

where a = , B(x,y) is the Beta function and i> is the effective magnetisation dependent 
spin velocity. A" and N' are velocity and hence field dependent pre-factors. These expressions 
for the susceptibilites are valid as long T < (H—H&). This is because above this temperature 
the hamiltonian in Eq.(ll) obtained by linearising the fermi surface around k' F is no longer 
valid. 

Though we have explicitly calculated the susceptibilites only in the IC regime for low 
temperatures, it is nonetheless interesting to study the behaviour outside this IC phase. For 
example, one can study how the system crosses over to the high temperature classical limit 
as a function of the parameters of the theory. We present a phase diagram as a function of 
T and H in Fig.3. This phase diagram is true only for H < H c2 because at H c2 the system 
makes a transition to another gapped phase. For temperatures smaller than the exchange J 
one obtains four regimes as indicated in Fig. 3. All the lines indicate crossovers and are not 
phase transitions. This phase diagram can be easily understood within the fermionic picture. 
Note that for the fermionised version of spin chains with a gap the complete dispersion for 
the excitations has the form Uk = ±v J 2 k 2 + A 2 . These two branches constitute the lower 
and upper bands with a gap 2A separating the two bands. With the full dispersion, the 
magnetisation is given by a function which depends on two parameters ^ and ~. The 
different regimes to be discussed below will be characterised by the varied behaviour of this 
function. For all T < H c \ — H one gets the gapped phase where the chemical potential i.e, H 
lies in the gap and the temperature is still small enough and does not excite particles to the 
upper band. This phase is characterised by an exponential decay of all correlation functions. 
When H > H c i and T < (H — H&), the chemical potential lies in the upper band and we 
get the the incommensurate or Luttinger liquid I regime as was already discussed in Sec. 
IIIB (cf. Fig.3). The dynamical exponent is z=l here. The susceptibilites in Eqs. (17-19) 
are valid in this Luttinger liquid regime. 

If the chemical potential lies at the bottom of the upper band and T ~ or equivalently 
when temperature is such that it can excite particles at the fermi level to the bottom of 
the upper band, one obtains the quantum critical regimei! defined bvlif — H&\ < T. In 
the quantum critical region the physics is governed by the fixed pointed at H = H& where 
K has a universal value -j. A simple way of understanding the behaviour in this region is 
presented below. Firstly, the curvature of the bottom of the band becomes very important 
in that the effective dispersion seen by the excitations is quadratic. Here one can expand 
u)k in A to obtain u® c = A + ^ in the quantum critical region. The dynamical exponent 
z = 2 in this region. A direct consequence of the quadratic dispersion is that 4 and It are 
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no longer independent scaling variables and only occur in the combination H and the 
magnetisation m which is just the number of fermions in the upper band has the scaling 
form m(T) = T^f{^^). Analo gous scaling functions exist for the spin-spin correlations 
also. The scaling functions for the correlators are functions of j, and H ~^ cl . Note that ^ 
and y are no longer independent scaling variables and only occur in the combination ~ cl . 
The spins are complicated functions of the fermions and this results in the scaling functions 
being difficult to obtain for reasons described in Ref. The quantum critical behaviour 
holds only as long as the quadratic form of the dispersion assumed is valid or in other words, 
only as long as temperature is such that the bulk of the excitations is confined to the bottom 
of the upper band. Similar phases were obtained in Ref. |23] in the context of spin-1 chains in 
a magnetic field, using a phenomenological theory of magnons with repulsive interactions. 
Here starting from a microscopic description, we find that these phases are generic to all 
gapped chains in a field and we have also obtained the exponents in the Luttinger liquid 
regime I. 

Now if temperature is increased further, there is a crossover to a fourth region. This 
regime was not obtained in Ref. |23] because the form uj^ was assumed for all values of k. 
Above a certain temperature excitations to higher k states occur and one starts probing the 
deviation from a quadratic dispersion and the k 2 approximation is no longer valid. As a 
result the quantum critical scaling for the magnetisation etc., will no longer be valid above 
these temperatures. The system then crosses over to a new region where temperature is 
large enough such that in addition to the excitations which now involve states way above 
the bottom of the upper band, there are transitions between the lower and the upper bands. 
This implies that the gap becomes irrelevant and the effective dispersion becomes linear 
in k. This crossover should occur for temperatures T of the order of twice the gap A but 
much lesser than the exchange couplings such that the interaction becomes irrelevant and 
the resulting behaviour is that of a Luttinger liquid. We denote this regime as the Luttinger 
liquid II. Classical behaviour sets in for temperatures greater than the exchange couplings. 
The widths of all these regions are dictated by the values of the gaps and the exchange 
coupling. For instance, the width of the quantum critical regime is fixed by the gap alone, 
and for magnetic fields away from H c \ this width is quite small and not as large as purported 
to be in Ref. ESI As a consequence, the scaling arguments apply only in the vicinity of the 
fixed point at H c \ where the width of the quantum critical regime is appreciable and not for 
fields far away from it. In contrast, the width of the Luttinger liquid II regime is given by 
the relative sizes of J and A. This width can be increased or decreased by tuning the ratio 
■j. Applying these arguments to the case of spin-1 chains where the gap of of the order of 
the exchange coupling J, one finds that for small fields the classical limit sets in soon and 
there are no sharp crossovers between the quantum critical, Luttinger liquid II and classical 
regimes. 

In the following paragraphs we discuss the physical significance of the susceptibilites 
calculated above. The susceptibilities are directly relevant to inelastic neutron scattering 
measurements where apart from certain magnetic form factors, the intensity is proportional 
to the (q,u) fourier transform of the full spin-spin correlator (S(x, t)-S(0, 0)). From Eqs.(17), 
we can see that the scattering intensity in the IC phase due to the S z correlator , obtained 
as a function of u for some fixed q, should in addition to the peak generated by the massless 
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excitations near Q = contains an extra peak corresponding to the massless modes at the 
incommensurate Q = 7r(l — 2m). Similarly, for the correlations perpendicular to the field, 
peaks are seen at the staggered mode Q = n and at Q = 2irm. These peaks are divergent 
at T = because of the power law correlations present at these values of Q. At finite 
temperatures the peak heights are finite and are determined by the T dependence of Eqs. (17- 
19). Such incommensurate features have been observed in inelastic neutron scattering data 
from copper benzoate which is a s = ^ system in a magnetic fiekfi Another probe is electron 
spin resonance (ESR) which can be used to directly probe the nature of the continuum of 
excitations at Q — 27rmil. 

One interesting question is whether the propagating modes corresponding to the uniform 
and staggered magnetisation are damped at finite temperature or in other words, is there 
any thermal broadening? This can be answered by studying the temperature dependence 
of the imaginary parts of the corresponding susceptibilites. First consider the case of the 
isotropic spin-| chain in zero magnetic field. Here it is known that the uniform magnetisation 
component of the S z correlator diverges as -\. The corresponding uniform and staggered 
susceptibilities at finite temperatures are given bj@il, Eqs.(17) with m = and K — |. 
Isotropy results in the same expressions for the perpendicular susceptibility. In this context 
v is the spinon velocity. The dependence of X®~ n on temperature tells us that the Q = it 
mode is damped at finite T. The temperature independence of the X®~° implies that there 
is no damping of the Q = mode at finite T@. This can also be understood from the fact 
that the exponent (3 in Eq.(18) is zero at K = \. There is no damping of this mode within 
the continuum approximation where the fermionic dispersion was assumed to be linear. 
Taking into account a small curvature of the fermion dispersion spectrum does not alter this 
result because the height of the peak which is given by the lifetime of the quasiparticles in 
the fermionic picture becomes infinite as q — > even at finite temperatures. Another way 
of understanding this absence of damping of the uniform mode is by noting the fact that 
the total magnetisation in all the directions commutes with the hamiltonian which in turn 
implies that there is no damping of this mode in the lattice spin system. 

We now discuss the damping of the various modes in the IC phase. The modes are the 
uniform magnetisation mode Q = 0, the staggered mode Q = it and the incommensurate 
modes at Q — l-nm and Q = 7r(l — 2m). From Eqs. (17-19) we see that the dominant 
contribution to the Q = n and Q = lum modes is from the perpendicular susceptibility. 
As in the Heisenberg case discussed above, the staggered mode at Q = tt has a damping 
factor proportional to T 2a ~~ 2 . Similarly the presence of a non-integral exponent (3 in x±~ 2nm 
results in the Q = 2irm mode being damped at finite temperatures by a factor proportional 
to T 2/3 . Similarly the behaviour of the Q = — 2m) mode is dictated by ^-j3- 7r ( 1_2m ) j n 

Eq.(17) and has a damping factor T 2K ~ 2 . We now consider the Q = mode. The dominant 
contribution to the damping of this mode arises from X^~°- From Eqs. (17), we find that 

the parallel susceptibility does not damp the Q = modes because X^~° is still independent 
of temperature. Because the system is not isotropic we need to check whether there is a 
sub-dominant contribution from the perpendicular susceptibility which damps this mode. 
To do this we first study the damping for other values of Q between and 27rm. Note that 
this damping can be studied at the resonance frequencies u = v(Q — q) or away from 
resonance. To obtain the leading temperature dependence of the damping of these modes 
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at resonance, we substitute q — 2-nm = Sq in Eq.(18). The resonance frequency to is fixed 
at v5q. The low temperature behaviour of the imaginary part of the susceptibility is given 
by 

^-lrpP+i + l^T?- 1 (20) 

This expression implies that there is a thermal broadening of the Q = or uniform magneti- 
sation mode i.e., 5q = —2ixm. Nevertheless, using the fact that the total magnetisation in 
the z— direction still commutes with the hamiltonian even in the presence of the magnetic 
field and that the time evolution of Y^i Sf involves only J2i $t we conclude that there is 
no thermal broadening of the uniform mode (Q = 0) mode. If u is not at the resonance 
frequency then for low T, there is no obvious thermal broadening and Imxj_~ 2nm is propor- 
tional to [v5q(v5q — 27rm)] /3 ~ 1 . This expression is not valid for uo — v5q < T and fails in 
the vicinity of 6q = 2irm. From these calculations we see that though bosonisation describes 
the physics correctly near Q = it and Q = 2-nm for the perpendicular correlations it does 
not describe modes far away from these two points well. This is not surprising in view of 
the fact that by linearising at k' F we take into account a lot of spurious states at the bottom 
of the band. To summarise, we find that the Q = mode is undamped whereas the gapless 
modes at Q = 2irm, ir(l — 2m) and tc are all damped. 



V. SPIN-i LADDER 



Another system which exhibits a behaviour akin to the systems studied above is the spin 
ladder. Here we consider a ladder with two identical and isotropic chains. Let the spins in 
chain 1 be labelled Si and that in chain 2. So. We bosonize this system in the same manner 
as above. We refer the reader to RefHJeZEa for details. We adopt the notations of these 
references and introduce the symmetric (triplet) and antisymmetric (singlet) combinations 
of the fields: c/) Sja = < ^ >1 ^ 2 and their respective duals 6 S and 9 a . The hamiltonian for the 
ladder is 

n = n a + n s (21) 

where 

H a = ^ J dx[Ku(7lU a ) 2 + (-|)(^0 a ) 2 + 9l COs( V / 80a) + ^2 COS(V26 a )} 

ft s =-L f dx[Ku(irTl s ) 2 + (^-)(d x <J) s ) 2 + g 3 cos(v / 80 s )] (22) 

Z7T J K 

where K — |, g\ — g% — g 3 — and J± is the interchain coupling and A some constant. All 
the cosine operators are relevant operators of dimension one. Eqs.(22) have the same form 
as the hamiltonian for an isotropic spin-1 system written in terms of two spin-| operators^. 
We note that the spin-1 system has fixed values for gi,g2 and g% and has no analog of a 
tunable parameter like J±. Nevertheless the results to be derived below apply to the spin-1 
system in a magnetic field. From Eqs.(21) and (22) we can infer the existence of gaps for all 
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non-zero values of J± in the spectra of both the fields <j) a and <j) s . The ground state of the 
ladder is a spin singlet and there exists a gap to the triplet excited state characterised by 
(j) s . Therefore, analogous to the dimerised chains ,we expect the vanishing of the gap and 
the onset of a gapless incommensurate phase for some critical value of the magnetic field. 
Since the magnetic field acts on both chains equally, we can easily see from the bosonization 
formulas that H affects only the <p s field 

n s ^n s + V2Hd x <f) s (23) 

Note that this hamiltonian has the same form as Eq.(lO). Therefore, using the results 
obtained in Sec. Ill, we can immediately see that the S field becomes massless while in the 
antisymmetric sector the 9 a field acquires a non-zero expectation value and the <fi a still has 
a gap. As a result the correlations of the <p a field decay exponentially. The value of K 
for the (p s field at H = H cl is K = ^. From Eq.(13),we find that K = | above H c \ also. 
This is because sinh(2#) = for the ladder. However, one should note that this is valid 
only in the region close to H c i and provided that the gaps and the magnetic field are small 
compared to the intrachain exchange coupling. Using the above results we can compute the 
spin correlation functions. There are two kinds of correlations: correlations within a chain 
and between chains. Here again the magnetic field changes the structure of the correlation 
functions in a manner analogous to that described in Eqs.(14). We summarise the results for 
the various correlators below. Since we are concerned only with the asymptotic behaviours 
we present only the power law contributions to these correlators. 

W (,, t )flJ(0,0)>=m» + g L[^ + ^] (24) 

Here r, t denote the chain labels. Unlike in the single chain case, in the ladder the alter- 
nating part at Q = 7T which is now shifted to Q = 7r(l — 2m) by the magnetic field decays 
exponentially. Similarly 

(S+ (x,t)Si (0,0)) = (S+(x,t)Sz(0,0)) 

= (-lf(x 2 -vh 2 yiK 
(Sf(x,t)Sz (0,0)) = t(-l) x (x 2 -vH 2 )-^ (25) 

Here again the uniform component is shifted to Q = 27cm and decays exponentially. The fact 
that (p a is massive results in an exponential decay of all the incommensurate contributions 
to the correlations. This exponential decay of all the incommensurate correlations in the 
ladder is in contrast to the single chain systems studied in the previous sections. Also note 
that except for certain exponentially decaying corrections, the interchain and intrachain 
correlators have essentially the same asymptotic behaviours. As already mentioned, the 
spin ladder in the magnetic field has the same exponents as that of the spin-1 chain in a 
magnetic field@§il. The dominant contribution to the perpendicular susceptibility for the 
ladder has the same form as that given in Eq.(19) with a = i and that to the parallel 

susceptibility by X\\~° °f Eq.(17). At the critical point H = H cl , the exponent a = \ for 
the ladder as well as the dimerised model. However, there is one big difference between the 
two models. In the dimerised model, the incommensurate parts of the (S Z S Z ) and (S + S~) 
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also show power law behaviours whereas, in the ladder they decay exponentially. This has 
a serious consequence for neutron scattering intensities. This is because for the dimerised 
system, at T = the power law divergences of the incommensurate parts of the dynamic 
correlations will result in a divergent peak at the incommensurate wave vector Q. On the 
other hand the exponential decay of the incommensurate correlations results in much smaller 
peaks at incommensurate Q whose finite height and width are determined by the gaps in 
the 4> a field. Though away from the critical point the exponents for the two models are no 
longer identical, the discussion presented above for the neutron scattering still holds. 

VI. NMR RELAXATION RATES 

With the help of the susceptibilites derived above, we can easily compute various quanti- 
ties that can be studied by neutron scattering and nuclear magnetic resonance (NMR). Here, 
we focus on NMR and in particular the the spin lattice relaxation time T\. The dominant 
contribution to T\ comes from the coupling of the nuclei to the lattice spins. Therefore, it 
is a good probe to study the nature of the lattice spin system. To obtain the temperature 
dependence of T\ we use the following formula in terms of the local susceptibility to calculate 
the sameS 

1 2ksT f dq 

— = h^o-p- / —FijiqjXijiq^iT) (26) 

Here the Fij are hyperfine form factors and Xij nas been defined in Eq.(16). In general 
these form factors are diagonal in i,j and do not vary much with q. For a system of non- 
interacting spins, (TiT) -1 is a constant. For interacting spin systems, the dependence on 
temperature could be more complicated because the underlying magnetic order plays a very 
important role in that it changes the effective magnetic field seen by the nuclei. Examples 
are the isotropic Heisenberg model where, T\~ l goes to a non-zero value as T — > and the 
spin-Peierls system where T{~ goes to zero at T = because of the gap to spin excitations. 

We now use the results obtained in the previous section to calculate the NMR rates. We 
first note that the magnetisation m only shifts the resonance frequency and does not change 
the form of the expressions for T-f ■ Depending on the kind of NMR done, one can probe 
specific correlations. This is especially useful for anisotropic spin systems and also isotropic 
systems in a magnetic field where the perpendicular and parallel susceptibilites are different. 
For instance, if the NMR was done on the nucleus of the lattice spin, then the relaxation 
occurs through a contact interaction and T-f 1 depends on x± alone. On the other hand if 
it is done on other neighbouring nuclei in the compound, the relaxation is through dipolar 
interactions and Tf 1 depends on x± and X\\- An amalgam of the two methods will be useful 
in isolating the two susceptibilites experimentally. 

Substituting the expressions for x± and x\\ derived in Eqs. (17), (18) and (19) in Eq.(26), 
we find that a straightforward power counting yields the following leading low temperature 
behaviour for the single chain models: 

[_]«njrJec/win _ A+T^' 1 + A||T 2 ^ -1 + B\\T (27) 
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where K is given by Eq.(13). Aj_,Am and Bn are constants independent of temperature. The 
suffixes _L and || refer to the contributions from the perpendicular and parallel susceptibilites. 
It is easy to see that the staggered susceptibilities dominate in Tf 1 . For the ladder model, 
Tf 1 is given by 

[ 1 | (odder = ^ ±T ^-1 + BT (28) 

The contribution coming from the A\\ term has not been explicitly written because it goes 
to zero exponentially as T — > 0. 

The temperature dependences of these rates for the dimerised and frustrated models at 
the critical point H ci are given below: 
Dimerised (5): 

±- = (A 5ll +B sll )T + A s± T^ (29) 

Frustrated models (J 2 ): 

±- = A nT -h + {Af± + B fll )T (30) 

In Table I, we present the leading low T contributions to Tf 1 at the transition, for the 
three models considered in this paper. There are other temperature dependent contributions 
to the NMR rate, but these go to zero at T = 0. As mentioned earlier there are two 
possible scenarios. One is that the nucleus probed does not correspond to the spins and 
the interactions are dipolar. Here at H = H cl , the Tf 1 diverges as T~2 for the three 
models. However, the divergent behaviour at low temperature in the frustrated model 
arises from the parallel susceptibility whereas in the ladder and dimerised systems it is the 
perpendicular susceptibility which leads to the divergent behaviour. This feature can be 
used to differentiate between the models as will be discussed in the following paragraph. 
As H is increased, we can see from Eq.(13) that K increases for the frustrated model and 
the divergence becomes weaker. Coincidentally K decreases for the dimerised system and 
the divergence of Tf 1 becomes weaker too. For the ladder, K does not change with H and 
the divergence persists and one has to go to higher fields to see a deviation from the T~2 
behaviour. If for some value of H, K decreases to | in the dimerised model and increases 
to | in the frustrated system , we see that Tf 1 does not diverge and Tf 1 goes to a constant 
A\\ +A± as T — > 0. This behaviour occurs because at some point the magnetic field becomes 
large enough such that the interaction which generates the gap becomes unimportant and 
should recover the exponents for the chain without the interaction i.e., the Heisenberg chain 
where Tf 1 is a constant. For fields greater than this value of H, the exponents vary like 
those of a Heisenberg chain in a magnetic field. The exponent approaches that of the XX 
chain or free fermions for sufficiently large fields. 

Another way of differentiating between the three models is if the NMR involves the nuclei 
of the relevant spins. Here, only the perpendicular local susceptibilites matter and a very 
interesting picture unfolds. At H = H cl , Tf 1 diverges as T _ 2 for the dimerised and ladder 
models whereas it approaches zero linearly in T for the frustrated model. Naively, for gapless 
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systems, we would have expected T{ 1 to diverge or go to a non-zero constant as T — ■> as 
it does in the case of the Heisenberg model ^1. This T dependence in the frustrated model is 
indeed strange because it is reminiscent of the T-f rates for spin chains with gaps where T-f 1 
also goes to zero as T — > but exponentially ! A divergence will be seen in the frustrated 
system if H is such that K > |. For increasing H the T~a divergence survives in the ladder 
but for dimerised systems Tf 1 becomes less and less divergent as T — ► and later saturates 
to a constant A± at T = for some value of H. This is very similar to what was seen in the 
case of dipolar interactions discussed above. However, note that the saturation values are 
different in both cases. We also note that the above discussion rests on the fact that these 
fields are smaller than which need not necessarily be the case. As mentioned above, 
for magnetic fields close to H C 2 where the ground state is nearly ferromagnetic, we expect 
the system to approach the free fermion limit i.e., K — > 1. As a result Tjf is expected to 
diverge as T~2 at H — for all the models irrespective of their values of n. 

Eventhough both the ladder and the dimerised systems have the same divergence at the 
critical point, one can differentiate between these two systems by studying the non-divergent 
contributions to T-f . For example, for the frustrated model, there exists a correction to 
Tf 1 proportional to T 2K+ 2k^ 1 . This exponent changes with increasing magnetic field and 
corrections of a similar nature do not exist in the ladder. These corrections should manifest 
themselves at not too low temperatures. However, inelastic neutron scattering should be 
able to differentiate between them as previously discussed in Sec. V. At temperatures large 
compared to the exchange couplings and the gap, from the analogy with fermions we expect 
that Tf = cste for all the models studied in this paper. A similar behaviour should be 
seen in the gapped phase i.e., H < H cl also. To summarise, we find that the kind of NMR 
experiments done can result in drastically different T-f for the three models. 

The results derived here can be checked in CuGeO^ and various other quasi ID systems. 
However, with most compounds being 3D the results obtained in this paper are applicable 
only in the temperature interval where the compound is effectively ID and that there is no 
3D magnetic ordering. Such a magnetic ordering in 3D could also result in divergent NMR 
rates. For instance, the onset of 3D Neel order at a certain temperature T/v, results in Tf 1 
diverging as (T — T/v) _ 5 for T > T/v0. This behaviour is valid for a temperature range of 
size T N . Since this divergence is the same as seen in the ladder, dimerised and frustrated 
systems at H = H c i it is important to establish whether such a divergence arises from the 
quasi-one dimensional or the 3D nature of the compound. This can be checked by working 
in the appropriate temperature interval where there is no 3D ordering or by increasing the 
magnetic field. If the compound is in the ID regime, its exponents vary with the field as 
predicted above and it if has 3D order the exponents do not vary with the field. 



VII. CONCLUSIONS 

We have studied the behaviours of various spin-| models in the gap less IC phase induced 
by an external magnetic field. For a generic gapped spin-| in a magnetic field, it was 
shown that the magnetisation is zero below H c i and rises as a square root above it. We 
found that the gapless behaviour in the IC regime is determined by the dimension of the 
cosine operator and hence different systems have drastically different properties. The results 
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presented here were obtained from a microscopic theory and not from a phenomenological 
theory as was done in the case of the spin-1 chain!!. We then discussed the implications of 
the finite magnetisation for the correlation functions. We found that the effect of the finite 
magnetisation was to shift the Q = it mode in the S Z S Z correlators to Q — 27rm and the 
Q = mode in the S + S~ correlators to Q = 7r(l — 2m). We also calculated the unequal 
time correlation functions and have provided explicit formulas for the various susceptibilites 
as a function of T,u,q and the magnetisation m. These were used to study the thermal 
broadening of the various modes in a single chain. We find that the modes at Q = ir, 
Q = 2Tim and Q = 7r(l — 2m) are broadened at finite temperatures whereas the Q = mode 
is not. Using the susceptibilites we also showed that neutron scattering intensities had extra 
peaks arising from the incommensurability in single chains but not in the ladder systems. 
We have also calculated the NMR relaxation rates as functions of temperature and have 
discussed the results for the dimerised,frustrated and ladder systems in detail. Using the 
fermion analogy, we find that the phase diagram for a generic gapped chain as a function of 
field and and temperature has five regions as shown in Fig. 3. In contrast to Ref. [23| where the 
system stays in the quantum critical regime for a wide range temperatures before it crosses 
over to the classical high temperature limit, we find that the system crosses over from the 
quantum critical regime to a second Luttinger liquid regime before it becomes classical. 
This intermediate temperature behaviours follow from the form of the dispersion spectrum 
of fermions with a gap due to interactions. Lastly, eventhough we have given the exact 
values of the exponents at the critical point H c \ alone, a knowledge of the magnetisation 
m from experiments can be used in conjunction with Eq.(13) to obtain the exponents and 
hence T± close to the transition at H&. We find that though magnetisation measurements 
cannot distinguish between the various models, techniques like NMR or neutron scattering 
which probe the dynamical spin-spin correlations will be able to do so. Other frequently used 
methods to study spin systems are EPR (electron paramagnetic resonance and Raman 
scattering. These methods might also be able to directly differentiate between the dimerised 
and ladder systems. We conclude by observing that it should be possible to verify the results 
obtained here in NMR measurements being done on CuGe0 3 and Cu2(C 5 Hi 2 N 2 )2ClJ^. 
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Table I 



Temperature dependence of T x at H = if cl 



Model 


K at H cl 


T -i 
J HI 


TTx 1 


Spin-Peierls (5) 


1 


T 


T-h 


Frustration (J2) 


1 

4 




T 


Ladder 


1 

2 


T 
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Figure Captions 



1. Schematic picture of the field dependent dispersion in the IC phase as seen by 
S Z S Z correlations (bold line) and S + S~ correlations (dashed line). Q = and 
Q = n are the usual commensurate modes and Q = 2nm and Q = n{\ — 2m) are 
the incommensurate modes. 

2. Plot of the magnetisation m (apart from an overall normalisation) versus magnetic 
field H for A = 0.4J at various temperatures T . H and T have been normalised 
by the exchange coupling J. The solid line corresponds to T = 0, the dashed line 
at T = 0.05 and the dashed dot line to T = 0.25. We see that even at sufficiently 
low temperatures the square root regime gets wiped out and m increases in a 
nearly linear fashion. 

3. Phase diagram of gapped spin-| chains in a magnetic field if as a function of 
temperature T. There are essentially 5 regions and all the lines indicate crossovers 
between these regions. 
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